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Abstract 

We present an application of the basic mathematical concept of complex 
functions as topological solitons, a most interesting area of research in physics. 
Such application of complex theory is virtually unknown outside the community 
of soliton researchers. 

Presentamos una aplicacion del concepto matemdtico de funciones complejas 
como solitones topologicos, una interesante area de investigation en fisica. Dicha 
aplicacion de la teoria compleja es prdcticamente desconocida fuera del circulo 
de investigation solitonica. 

Introduction 

The complex variable z = x + iy, where x, y £ 3ft and % = is one of the 

most familiar and useful concepts in mathematics, with a very large number of 
well-documented applications in science. 

Over the past few years some interesting nonlinear models in physics have 
received a lot of attention, models bearing the so-called solitons or energy 'lumps'. 
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Some of these models exemplify yet another important application of complex 
functions, with functions as simple as f(z) = z describing a soliton configuration. 
Unfortunately, despite the vast literature dealing with complex analysis plus 
applications, one finds no mention of the starring role of z as a soliton. Reference 
to such an extraordinary role is found only in highly specialised research books 
and journals, hence the existence of z as a soliton field is practically unknown 
outside the group of specialists in the area. 

Using the nonlinear sigma 0(3) (or CP 1 ) model in two spatial dimensions, 
the present work illustrates the context in which z stands for a lump of energy. 
This fact is most remarkable and, given the growing importance of solitons in 
physics, we believe that more physicists should know about it. They will find 
this fresh utility of complex variables quite appealing. 



Complex theory is a very important branch of mathematics. As a brush-up 
we just recall that many integrals given in real form are easily evaluated by 
relating them to complex integrals and using the powerful method of contour 
integration based on Cauchy's theorem. In fact, the basis of transform calculus is 
the integration of functions of a complex variable. And intersections between lines 
and circles, parallel or orthogonal lines, tangents, and the like usually become 
quite simple when expressed in complex form. 

Familiarity with the complex numbers starts early, when at high school the ba- 
sics of z are taught. Then in college algebra/calculus one learns some more about 
complex variables, with immediate applications to problems in both physics and 
engineering like electric circuits and mechanical vibrating systems. Later on 
complex holomorphic (analytic) functions are introduced, and then applied to a 
variety of problems: heat flow, fluid dynamics, electrostatics and magnetostatics, 
to name but few. 

The concept of analyticity is extremely important. Many physical quantities 
are represented by functions f(x,y), g(x,y) connected by the relations d x f = 
d y g, d y f = —d x g, where d x f = df/dx, etc. It turns out that / and g may be 
considered as the real and imaginary parts of a holomorphic function h of the 
complex variable z: 



The equations linking / and g are the Cauchy-Riemann conditions for h(z) being 
holomorphic, and can be written compactly as 



When h is a function of z = x — iy, the complex conjugate of z, the condition 
(j2J reads d x h = id y h, and h(z) is said to be anti-holomorphic pQ. 

We hereby show how functions of the type Q describe solitons, giving yet 
another fundamental, if little known, application of analytic complex functions. 



1 Complex theory 




(i) 



id v h. 



(2) 
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2 Solitons 



Nonlinear science has developed strongly over the past 40 years, touching upon 
every discipline in both the natural and social sciences. Nonlinear systems appear 
in mathematics, physics, chemistry, biology, astronomy, metereology, engineer- 
ing, economics and many more [21 IS] • 

Within the nonlinear phenomena we find the concept of 'soliton'. It has got 
some working definitions, all amounting to the following picture: a travelling 
wave of semi-permanent lump-like form. A soliton is a non-singular solution of a 
non-linear field equation whose energy density has the form of a lump localised 
in space. Although solitons arise from nonlinear wave-like equations, they have 
properties that resemble those of a particle, hence the suffix on to covey a cor- 
puscular picture to the solitary wave. Solitons exist as density waves in spiral 
galaxies, as lumps in the ocean, in plasmas, molecular systems, protein dynamics, 
laser pulses propagating in solids, liquid crystals, elementary particles, nuclear 
physics... 

According to whether the solitonic field equations can be solved or not, soli- 
tons are said to be integrable or nonintegrable. Given the limitations to anal- 
itycally handle nonlinear equations, it is not surprising that integrable solitons 
are generally found only in one dimension. The dynamics of integrable solitons 
is quite restricted; they usually move undistorted in shape and, in the event of a 
collision, they scatter off undergoing merely a phase shift. 

In higher dimensions the dynamics of solitons is far richer, but now we are in 
the realm of nonintegrable models. In this case analytical solutions are practi- 
cally restricted to static configurations and Lorentz transfomations thereof. (The 
time evolution being studied via numerical simulations and other approximation 
techniques.) A trait of nonintegrable solitons is that they carry a conserved 
quantity of topological nature, the topological charge -hence the designation 
topological solitons. Entities of this kind exhibit interesting stability and scatter- 
ing processes, including soliton annihilation which can occur when lumps with 
opposite topological charges (one positive, one negative) collide. For areas like 
nuclear /particle physics such dynamics is of great relevance. 

Using the simplest model available, below we illustrate the emergence of topo- 
logical solitons and their representation as complex functions. 

3 The planar 0(3) model 

Models in two dimensions have a wide range of applications. In physics they are 
used in topics that include Heisenberg ferromagnets, the quantum Hall effect, su- 
perconductivity, nematic crystals, topological fluids, vortices and solitons. Some 
of these models also appear as low dimensional analogues of forefront non-abelian 
gauge field particle theories in higher dimensions, an example being the Skyrme 
model of nuclear physics [U [3] . 
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One of the simplest such systems is the 0(3) or CP 1 sigma model in (2+1) 
dimensions (2 space, 1 time). It involves three real scalar fields <pj (j=l,2,3) 
functions of the space-time coordinates (t,x,y) [B1E]- The model is defined by 
the Lagrangian density 

^ = iE[(^) 2 -(^) 2 -(^') 2 ] (3) 

where the fields, compactly written as the vector in field space <f) = (<pi, fa, fa), 
are constraint to lie on the unit sphere: 

S^ = {fa.^ = l}. (4) 

The Euler-Lagrange field equation derived from ©-© has no known analytical 
solutions except for the static case, which equation reads 

vV-(M)l=o [v 2 = d 2 x + d 2 y }. (5) 

The CP 1 solitons are non-singular solutions of (JHj). Without the constraint 
(J3J the said equation would reduce to V 2 (j) = 0, whose only non-singular solutions 
are constants. The condition (0} leads to the second term in equation which 
does yield non-trivial non-singular solutions as we will later see. 

Solitons must also be finite-energy configurations. From @ we readily get 
the static energy 

3 

72, 



E = /^E[(^) 2 + (^) 2 ]rf 2 

3 

= /^E(v^)(v^)d 2 x [v=(d x ,d y )\ 

= — J (X7(f)).(X7(f>) rdrdO (in polar coordinates r, 8). (6) 

We ensure the finiteness of E by taking the boundary condition 

lim (j)(r,6) — ► 0o (a constant unit vector independent of 0), (7) 

since the integrand in © will thus tend to zero at spatial infinity : 



Urn r| V0| = Jhn r ^/ (<9 r <^) 2 + {^d e $) 2 -> 0. (8) 



3.1 The complex plane 



We are thus considering the model in the x—y plane with a point at infinity, ie, the 

(x) 

extended complex plane which is topologically equivalent to the two-sphere Sv; • 
The finite energy configurations are therefore fields (j> defined on ^Ujoo} = 
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and taking values on Sv, . In other words, our finite-energy fields are harmonic 
maps of the form -► S { ^ ] 0. 

(x) 

We may imagine the coordinate space Sv; as made of rubber and the field 
space 5g as made of marble; the map (f constrains the rubber to lie on the 
marble. Then with each point (x, y) in the rubber we have a quantity 

f = V 2 <£- (<£.VV)0 

at the point in marble representing the tension in the rubber at that point. 
Thus the map is harmonic if and only if <p constrains the rubber to lie on the 
marble in a position of elastic equilibrium, r = 0, which is just equation Ipj). 
These are our finite-energy configurations, of which the soliton solutions are a 
subset. 



3.2 Topological charge 



In general, as the coordinate z = (x, y) ranges over the sphere once, the 
coordinate (f = (</>i,02 5 03) ranges over N times. This winding number is 

(x) 

called the topological charge in soliton parlance, and classifies the maps — > 
into sectors (homotopy classes); maps within one sector are equivalent in 
that they can be obtained from each other by continuos transformations. 

An expression for the topological charge is obtained by expanding the coor- 
dinates 4>j of the area element of m terms of coordinates (x,y) in s!f \ and 
integrating off. In plainer language, from the college formula that computes the 
flux of a vector A through a region D of a surface S : 

/ A.iidS [n a normal unit vector], (9) 
Jd 

the topological charge N follows by calculating the flux of A = -^-j<f through 
the sphere of radius a = 1: 

A.ndS ^ [ —^$.$dS = N. 

D Js ( 2 x) 4vra 2 

Notably, a field with topological charge N describes precisely a system of N 
solitons. 

In order to actually find charge-iV finite-energy solutions, it is convenient to 
express the model in terms of one independent complex field, W, related to (p 
via the stereographic projection 

W = ^±^. (10) 

1-03 

In this formulation, the topological charge is given by 
„ If \d z W\ 2 - \d- z W\ 2 ^ 2 
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connected with the energy © through 

E>2ir\N\. (12) 

4 Lumps 

The solitonic solutions we seek correspond to the equality in (fT2l 9,|TO1^2. That 
is, in a given topological sector N the static solitons of the planar CP 1 model 
are the configurations whose energy E is an absolute minimum. Combining (|llj) 
with E = 2k\N\ we find that solutions carrying positive or negative topological 
charge satisfy, respectively, 

d- z W = -> d x W = -id y W, (13) 

d z W = -> d^W = id y W. (14) 

But recalling equation (J2J) we immediately recognise the above equations as the 
Cauchy-Riemann conditions for W being a holomorphic function of z or z. This 
is most remarkable. 

For instance, a single-soliton solution (N = 1) may be described by 

W = z [note that this satisfies equation (|13j)]: (15) 

its energy density distribution is given by 

1 + \z\ z 

Plots of (jl()|) reveal a lump of energy localised in space, as shown in figure^ The 
same energy corresponds to W = z, which has N < and sometimes is referred 
to as an anti-soliton. 

A more general N = 1 solution is given by a rational function 
W = X(z — a)/(z — b), which we should note is non-singular: W(z = b) = oo 
corresponds to = 1, the north pole of according to (|1U|) . A prototype 
solution for arbitrary N > is X(z — a) N . 

The dynamics of these structures is studied by numerically evolving the full 
time-dependent equation derived from ©, with the fields W(z) as initial condi- 
tions pa mi!- 

Sigma CP 1 -type models have several applications, noteworthy among them 
being the Skyrme model in (3+1) dimensions where the topological solitons stand 
for ground states of light nuclei, with the topological charge representing the 
baryon number. 

The role of complex functions as topological solitons deserves widespread 
attention and should not be missing from the modern literature dealing with 
complex theory and its applications. 
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Figure 1: The energy distribution corresponding to the soliton W = z. 
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